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Abstract 



o 

V^ It has been demonstrated that Lattice Bohzmann schemes (LBSs) are very efficient for Computational AeroAcoustics 

(CAA). hi order to handle the issue of absorbing acoustic boundary conditions for LBS, three kinds of damping terms 
are proposed and added into the right hand sides of the governing equations of LBS. From the classical theory, these 
terms play an important role to absorb and minimize the acoustic wave reflections from computational boundaries. 
P^ Meanwhile, the corresponding macroscopic equations with the damping terms are recovered for analyzing the macro- 

scopic behaviors of the these damping terms and determining the critical absorbing strength. Further, in order to 
detect the dissipation and dispersion behaviors , the linearized LBS with the damping terms is derived and analyzed. 
^ The dispersive and dissipative properties are explored in the wave-number spaces via the Von Neumann analysis. 

The related damping strength critical values and the optimal absorbing term are addressed. Finally, some benchmark 
problems are implemented to assess the theoretical results. 
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o 

f^ 1. Introduction 

>-. 

r^ In the past decades, there has been an increasing highlighted interest in developing and applying the LBS for 

Qh aeroacoustic applications [1, 2, 3, 4, 5, 6], because the LBS has been developed into an innovative mescoscopic 

numerical method for the computational modeling of a wide variety of complex fluid flows and acoustics [1, 2, 3, 
,—1 4, 5, 6, 7, 8, 9, 10]. In the field of the aeroacoustic applications, it has been demonstrated that the LBS possesses 

^ the low dispersion and low dissipation properties for capturing the weak acoustic pressure fluctuations [3, 4, 5, 6]. 

^^ Recently, the MRT-LBM (multi-relaxation time lattice Boltzmann method [8, 11]) has been improved for the acoustic 

^V_- application with the nice dissipation and dispersion relations. The improvement of the MRT-LBM overcomes the 

vQ drawback of the large bulk dissipation [ 1 ] . The improvement has made the MRT-LBM appearing as a very well-suited 

method for acoustic simulations. 
^;-J In the simulations of the realistic flows, artificial computational boundaries are defined around the flow region 

p^i of physical interest. The space outside the computational domain is neglected. In some cases, the region of interest 

T— I should extend to infinity, especially for the computational aeroacoustics [12]. Because of the use of a finite compu- 

L| tational domain and artificial boundary conditions (e.g. Dirichlet, Neumann or Robin types) the outcoming acoustic 

• ^H waves may be reflected. The reflected waves often have a significant influence on the flow field and may overwhelm 

^% the physical acoustic waves [12, 13, 14]. The obvious ways to suppress the reflected waves are to introduce artificial 

dissipation (by upwinding) or to increase the value of physical viscosity (or add hyperviscosity) [12]. The simplest 
way is to define sponge/absorbing layers is to add an additional penalty term to the governing equations to suppress 
the computed solution and to match a prescribed or precomputed solution [13, 14, 15]. The most used and efficient 
forcing term is given by [14] 

-/r(x)(?-?i-ef), (1) 
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where q denotes the physical variables under consideration. Physically, the term (1) plays the role of a linear friction. 
Generally, the parameter ;t' , called the sponge strength, is a function that depends on space only and ^ret is a function 
of space and time (especially, ^lef denotes any variable of the far fields or time-dependent mean value). In most cases 
X will vary smoothly between zero in the physical domain and a positive value in the absorbing laryer [13, 14, 15]. 
Provided ;\f is large enough, in the classical theory of absorbing layers, the outgoing disburbances are exponentially 
attenuated when crossing the layer[12, 14, 15, 16]. The numerical solutions inside absorbing layers need not be 
physical as long as the use of the region does not introduce significant reflection back into the physical domain. 
Meanwhile, the absorbing layer is numerically stable. Compared with perfectly matched layers, the absorbing layer 
techniques can be applied to a wider class of problems and are often coupled with the characteristic boundary condition 
[16]. 

The definition absorbing layers for LBS is an emerging resarch topic. It must be emphasized that there is a clear 
need for theoretical researches on that topic, to improce the potential of LBS for CAA applications. From the macro- 
scopic equations standpoint, we are interested in the macroscopic variables which are obtained by the mesoscopic 
distribution statistics. We propose in the present paper several possible absorbing terms which are enforced on LBS. 
By the classical Chapman-Enskog multi-scale expansion, the macroscopic systems are recovered in the absorbing 
layer. With the aid of the macroscopic systems, the related analysis is implemented in spectral space. The emphasis 
is put on the coupling between cr(x) and LBS relaxation parameters, especially on finding the critical value for cr(x). 
The critical expressions of cr{x) will offer us a guide for choosing the reasonable absorbing strength. This is different 
from the classical absorbing layer theory in Navier-Stokes equations (NSEs). In the classical absorbing layer theory, 
cr can be regarded as a penalized parameters, and theoretically, cr can be chosen as any positive real number However, 
in order to avoid a stiff problem, cr will not be chosen very large. In LBS, cr is a finite parameter dependent on the 
relaxation parameters, as shown hereafter. 

In next section, the fundamentals of LBS are reviewed and the possible absorbing terms are proposed. In the 
third section, in the spectral space, the dispersion and dissipation relations are analyzed with the aid of Von Neumann 
method considering monochromatic wave solutions. In the fourth section, the numerical simulations are performed to 
validate the proposed absorbing terms considering some benchmark problems. 

2. Fundamentals of LBS and absorbing terms 

In this section, the fundamental theory of LBS is briefly reviewed. Then, we give the linearized LBM problem 
coupled with absorbing terms and establish the relation between the classical BGK-LBM and the Navier-Stokes 
equations in absorbing layers. 

2.1. Lattice Boltzmann schemes 

The governing equations of the lattice Boltzmann schemes are described by the following universal form [8, 11] 

fiix + ViSt, t + 5t)^ Mx, t) + Aij [ff'hx, t) - fj(x, t)), < /, ;■ < A^, (2) 

where v, belongs to the discrete velocity set 'V, fi(x, t) is the discrete single particle distribution function corresponding 
to v, and ff denotes the discrete single particle equilibrium distribution function. Generally, ff can be expressed 
by a combination of a linear part / ' '''' (x, t) and a quadratic part / '^'' (x, t) [9] 

ff'^(x,t)^t'"^(x,t) + ff-^'^(x,t). (3) 

6t denotes the time step and A^ + 1 is the number of discrete velocities. A,y is the generalized relation matrix. From 
here on, the repeated index indicates that the Einstein summation is used except for some special explanations. Let 
X.eW' id denotes the spatial dimension) denotes the lattice system, and the following condition is required [2] 

X + Vj5t € £, (4) 

that is to say, if x is a node of the lattice, x + Vjdt is necessarily another node of the lattice. Generally, for BGK-LBM, 
the relaxation matrix is given by 



A;,/ = s6ij, (5) 

where s denotes the relaxation frequency of BGK-LBM. If the relaxation matrix A is defined by 

A^M-^SM, (6) 

where 5 is a diagonal matrix which denotes the relaxation parameters of MRT-LBM, which is given by 

S ^dmg({0,0,0, Sd^i,...,S!, }). (7) 

d+l N-d-l 

M - (Mij) is the transformation matrix (see Appendix A for details of D2Q9 and D3Q15), which satisfies 

the following basic conditions [8] 

Moj^l,Maj^v'],(l<a<d). (8) 

The macroscopic quantities are defined by [2, 8] 

nii^M^jfj, mf''^^M,/p\ (9) 

By the simple algebra, the standard isothermal MRT-LBM is recovered in the following form [2, 8] 

mi^Wi^inf°'\0<i<d, (10) 

and 

mi{x + 6tVj, t + 6t)- nijix, t) + Sj (mf'^\x, t) - ot,(x, t)\,d+\ <i<N. (11) 

It is necessary to point out that for isothermal flows, the number of the conservative quantities is equal to li + 1. 
According to the work of Lallemand and Luo [8], the relaxation parameters in Eq. (11) should satisfy the following 
stability constraints 

Si&{Q,2), d+\<i<N. (12) 

2.2. Absorbing terms 

Following the theory of Israeli & Orszag [14], Eq. (2) coupled with the absorbing terms has the following form 

/;(x + Vi6t, t + 6t)^ /i(x, t) + A,v [ff''\x, t) - fj(x, t)) + Str^j {ff'\x, t) - f*(x, t)) , (13) 

where F^ is the generahzed absorbing strength defined by 

F = M-^I.M, (14) 

and /j '^ denotes the reference state of /) and /* denotes the possible representations of mesoscopic distribution 
functions. From the classical theory, a natural choice of/* is fj. The matrix S in the expression (14) is defined by 

S = diag{cro, . . . , cr;v), (15) 

where cri is the absorbing coefficient for each first-order moment of /. In the expression (14), the transformation 
matrix M is given in Appendix A. The choice of the generalized absorbing strength expression (14) is based on the 
ideal for the different modes, the different absorbing strength can be applied as indicated in [13]. For simplicity, we 
consider (Tq = . . . = ctj^ - x, and the matrix F become the following diagonal form 

r,7 = &ijx- (16) 

In this paper, the researches are focused on the absorbing strength with the form (16) . 



Generally, the reference distribution function / (x, f) is expressed by the time-averaged mean value [16]. Here, 
we define the expression of /."^ '(x, t) by the equilibrium distribution function as follows 

/j'-^f>(x, t) = ./^^'^ V^"''(x, 0, uf^'^f^Cx, 0,0, (17) 

where p^'^^^\x, t) and u'"'^'(x, t) denote the reference (or base) density and velocity, respectively. Especially, let 
p*"'^^(x, t) and u*'^'^'*(x, t) denote the far fields. The influence of the absorbing terms on the macroscopic equations 
can be interpreted as the source terms of the macroscopic equations. Based on this consideration, the corresponding 
source terms can be obtained by the classical Chapman-Enskog procedure. 

It is necessary to discuss the choices of / and /*. Intuitively, let / "^ and /* expressed by the equilibrium 
distribution functions. This consideration makes us easy to analysis the macroscopic behaviors of the absorbing 
layers. Meanwhile, from these choices, we only need to handle the macroscopic statistical quantities and reference 
quantities. Compared with handling the distribution functions, these choices are effective and simple. Especially, 
provided you consider the time-dependent mean quantities as the reference state, the time-dependent mean of the 
macroscopic statistical quantities is easier to be handled than the mesoscopic distribution functions from the view of 
saving memory and CPU cost. Another more careful consideration is that generally, for aeroacoustic problems, if the 
background flows or far fields are used, it is diflicult to construct a suitable / . If/ is defined by the equilibrium 
distribution functions, this will lead to some potential instabilities when /* is defined by /,. The analysis will be given 
in Sec. 3. From above considerations, there exist three-kind simplest absorbing terms. 



2.2.7. Type I absorbing term 

( 
ff^'ipi', uf, and /;(x, t) = fj(x, t) we have 



Combining Eq. (13) and Eq. (17) and using the far field as the reference state, and considering / (x, — 



Mx + Vi5U t + St)^ /(x, + s' {f''ip\ u\ - fiix, O) + dtFiipf, u-^,p% u*, 0, (18) 

where s' = i -h 6tx, Fj(p-^, v/ ,p*, u*, is defined by 

Fiipf, uf,p\n\ t)^x [ff''\p*\ u-^, t) - ff°'\p\v^\ t)) . (19) 

In Eq. (18), p* and u* denote the density and velocity, which are defined by the following parameterized forms [17] 

p* =2_^fi+ nSt 2_^ Fi, p*ul = 2_j ^iafi + "^<^f /_j CiaFj. (20) 

/ / / / 

Then, we have the following macroscopic equations (see Appendix B for details) 

d,p* + daip'uj = (1 H- ns')x(p^ - p*) + (i^dtipf - p*) + i^,„xda {f^u{ - p*ul^ , 
. dtip'ul) - d/}(v'p*(dau*fj + dpuj) + dp [p*ulup + p*6afs) = (1 + ms')x [p^ui - p*u^ + (21) 

(mXdi {p^ui - p*ul) - (o- + -:r\xdp [p^ uWp + P^ Sap - p*ulup - p*5ap) , 
where s' , v', cr, ^„ and ^,n are defined by 

s' ^s + 6tx,V ^clij,- y '^f''^ = (j, - ^J-J^^" = «^M 1 - y - f ' ^m = m5tU - ^ - f ■ (22) 

According to the definition of the effective relaxation frequency s' , it is known that for Eq. (18), an effective 
negative viscosity v' will appear, \f x satisfies the following condition 

s' ^ s + 5tx> 2. (23) 



2.2.2. Type II absorbing term 

According to the deviation in Sec. 2.2.1 and considering /''''^(x, i) = /■'"'^{f/, u-^, t) and /*(x, t) = f^^\p\ u\ t), 
we propse the second-type absorbing term for LBS 

/■(x + Vi6t, t + 6t)^ /(x, + s {f''\p*,n\ t) - /(x, t)) + dtPf^Xpf , uf ,p\u\ t), (24) 

where Fiif/ , u^,p*, u*, f) is defined by 

Ff\f^\ nf,p\n\ t)^x U'^'^V, u^, " ff'^P, u*, O) , (25) 

In Eq. (24), p* and u* denote the density and velocity, which are defined by the following parameterized forms [17] 

p*=^/; + «ft^Ff^\ pul^Yj'=i,fi+m5tY,CiaFf''\ (26) 

The corresponding macroscopic equation is similar to Eq. (21). It reads (see Appendix B for details ) 

'' d,p* + da{p*ul) = (1 + ns)x(p^ - p*) + CnXdtiP* - P*) + imX^a [p* ui - p* U^) , 

d,(p*ul) - d/}(vp*(dau*fi + dpuj) + dp [p*ulup + p*5ai3) = (1 + ms)x [p^ui - p*ulj + (27) 

1 {^«, (M - P-<) - (.. . f ).% (M4 . /■>„, - P-<«i - P-s.„) . 

where cr, („ and ^„, are defined by 



■=4'^-\h^-[l-\h^"-'''i'-^)-T^'"-'"'i'-'^-T 



(28) 



The diff'erence between Eq. (18) and Eq. (24) lies in the relaxation frequencies s' and s. In Eq. (24), the effective 
relaxation frequency s is kept in the original form while, in Eq. (18), the original relaxation was modified, possibily 
leading to a negative viscosity if 6tx is too large. This problem is now cured. 

2.2.3. Type III absorbing term 

In the recovered macroscopic equations (21) and (27), in the left hand side, there exist momentum-velocity coupled 
terms p*M>^ in recovered equations. Considering /j'"'^(x,r) = f'j^''''^\f/,u-^,t) and /*(x,f) = /|^''"'V%u*,0, these 
terms can be dropped by the following third type absorbing term 

fi(x + vi6t, t + 6t)^ fi(x, t) + s [f^hp^ix, t), u*(x, 0, t) - fi(x, t)) + 6tF^^- ^^' {pUx, t), u^ix, t),p*{x, t), u*(x, t), t)) , (29) 

where u*(x, t) is the equilibrium velocity. F ' '''' (u/(x, f), u*(x, f), f) is defined by 

f'-I^' "'\pf, nUx, t),p\n*(x, t), t)^x [f?^' ''V^, u^(-^, 0, - t' ''V% u*(x, t), t)) . (30) 

In Eq. (29), p* and u* denote the density and velocity, which are defined by the following parameterized forms [17] 

p*^Y. f' + '"^^ Z ^^ ''*' p*< = Z '-'-f- + ""^^ Z ^-^^ '"■ ^31^ 

Then, the following macroscopic equations are obtained (see Appendix B for details) 

i d,p* +daip*u*a) = (1 + ns)x(p^ - p*) + ^„xd,(p^ - p*) + ^„{>(da[p^ u^ - p*ul) , 

\ d,(p*K) - d/}(vp*(daU*^ + dpuD) + dfi {p*ulu*^ + p*6ais) = (1 + ms')x [p^ui - p*ulj + ^mX^t [p^ui - p*ul] , 

(32) 



where the parameters v, ^„ and (,„ are defined as follows 

'1 1 



St, („ = ndt 



('-I) 



6t 

2"' 



im = m6t 



('-f) 



ft 

2' 



(33) 



In this section, the fundamental theories of LBS are reviewed. Based on the classical absorbing layer strategies, 
the three types of absorbing terms are proposed for LBS, and the corresponding macroscopic equations are given. In 
the following section, the behaviors of the absorbing terms are studied. 

2.3. Profiles of absorbing strength cr 

In the absorbing layers, the profile of cr could not be specified as a uniformly distributed parameter. A uniform 
o" distribution will induce a significant wave reflection from the interfaces between the wave propagation domain and 
the sponge domain [13, 14]. The most popular profiles of o" are listed in Table 1. In the current research, the first kind 
of cr will be used which was first proposed in [14]. For this kind of absorbing profile, cr{x) satisfies cr(xo) - o-(L) = 0. 
In the original paper of Israeli and Orszag [14], by taking A - n - A, they obtained the best results for the acoustic 
wave equation. In order to handle the current problem, we given following normalized cr^ for A = n = 4 as follows 



_ 3125(L-x)(jc-xor 
^^~ 256(L - xo)5 
The normalized a-(x) profile for xq -Q and L = 1 is given in Fig. 1 



(34) 




Figure 1: The normalized (t(x) profile: A = n = A, xq = 0, L = \. 



3. Analysis of LBS coupled with absorbing terms 

In this section, the influences of the absorbing terms on the dissipation/dispersion properties are explored looking 
at the linearized forms of the recovered macroscopic equations and LBS. 

3.1. The linearized macroscopic equation with the absorbing terms 
If we consider n - m, Eqs. (21) and (27) can be expressed by 



dip* + daip'uj = (1 + ns)x(p^ - p") + ^nXdtip^ - P*) + (.,{>(da {f^u{ - p*ul) , 
d,{p*ul) - di}{vp*{daUp + dpuj) + dpipulup + p*6a/s) = (1 + ns)x[p'ui -p''ul) + 

^nXdt {p^ui - p*ul) - (o- + y j.r5/J {p^uWfi + P^Sa/i - P* U^Uf^ - p*6alj) , 



(35) 



Table 1: Newtonian cooling function: possible profiles of o" from the literature [14, 18, 19]. 



Damping Types 



cr(x) 



Parameters 



Polynomial 



Reyleigh 



(x - xq)"{L - x){n + !)(« + 2) 



(L - xq) 



n+2 



Tl- 



X- Xq 



[L- xq) 
a[l - cos(87r(x - xq - 1))] 



xq: Position of sponge layer; A: Absorb- 
ing strength; L: Width of sponge layer; for 
most choices, n - A. 

Tj: Absorbing strength; L: Width of sponge 
layer; n - 3. 

Xq: Position of sponge layer; a: Absorbing 
strength. 



where v denotes v or v', s denotes s or s', and Eq. (32) can be rewritten as follows 



d,p* + da(p*ul) = (1 + ns)x(p^ - P*) + ^nXdtiP^ - P*) + (nX^a [p^ »{ - p* U^j , 

dtip*iO - df}{vp*{daUp + dpuj) + dp [p*ulup + /?*(5„/j) = (1 + n% [p^u{ - P*u^ + inX^t {f^u{ - p*u^ , 



(36) 



In order to detect the dissipative behavior, we consider the following simplified ID case of Eq. (35) (the flow 
propagation is only along x-axis) 



d,p* + dAp'ul) = (1 + nS)x(p^ - p*) + t^nX^tif^ - P*) + l^nX^x {p'ui - P*u^ , 

d,{p*u\) - dAvp*{dxul + d_yj) + d^ (p*u*X-c + P*) = (1 + «% {p^u{ - p*u^ + ^„xd, {p^u{ - p*u^ (37) 

-{0-+ 2Jxdx {p^uiu{ + pf - p*uy^ - p*) . 



We consider the foUowing decomposition of p* and m* 



P* =P° +P\Ux = "? + "x 



(38) 



For the sake of simplicity, we take u{ - u^ - andp" - pf - const - 1, leading to the following Unearized system 
d,p' + dy, = -(1 + ns)xp' - ^nXdtP - (nX^xU^, 
Oy - Ivdlu'^ + c\d^p' = -(1 + «%«', - i;,„X^tu'^ + c^ (cr + —\xdxP'- 

From the above equations, we obtain 



(39) 



A ■ ( 1 + UW,p + (Bx + A-Bx + B^„x')d,p - 2Avdi(d,p) - (2vBx + Acj - Ac\Cx)dip + BYp = 0, (40) 

where A - I + (nX^ B -\ +ns and C - cr + dtjl. 

Now, considering a monochromatic wave solution p' - pgexp[i(fef ■ x — (jDt)\, we get the following k- u) relation 

A ■ (1 + C„xW +i-(Bx + 2Avkl +A-BX + B(„x^)aj - IvBxkj - Ac^, - bY + Ac]Cx^, = 0, (41) 

Considering a monochromatic solution, the similar k - co relation for the linearized form of Eq. (39) can be 
obtained as follows 



A ■ (1 + i„x)^' +i-{Bx + 2Avki +A-BX + Bi„x')co - 2vBxki - AcX ' B'x = 0. (42) 

7 



If y vanishes in Eq. (39), then, we are left with the following ID form 

dtp' + 5^m', = -(1 + ns)xp' - t^nX^tP' - ^nXdxU'x, 

dtu'x + c]d^p' = -(1 + ns)xu',. - (,„xd,u'^ + c]\a- + -^jX^xP'- 

The above equation is different from the classical ID model 

d,p' + dXx = -CT-p(x)p', 
d,u'^ + c^dxP - -cr,i{iC)u'j^. 



(43) 



(44) 



Mathematically, the above classical absorbing strategy is similar to the penalty method [13]. As indicated, if 
(p', u'^) is suitably smooth and the norm ||(crp, cr„)|| suitably large, the spatial derivatives of (p' , u\) will be negligible 
compared with the right hand side. Eq.(44) and the solution decays exponentially with the time scale 7~ [13]. If f > 7~, 
the right hand side terms will not dominate Eq. (44) and the convergence property will depend on the initial condition 
and the coercivity of the convection operator. From the viewpoint of penalty methods, cTp and cr„ could be chosen as 
any positive large numbers. Now, let us rewrite Eq. (43) as follows 



d,p' + 5.fM', 



(1 +ns)x , 
' l+inX^ 



, {cl-cU(r + 6t/2)x) , (l+«% , 

a,Mj H ; ; OxP : U 



(45) 



We now consider h 
inequality is satisfied 



1/2 and 6t 



1 + inX 1 + ^nX 

\, („ - -5/4. It is clear that („ is negative. That means if the following 



the effective coefficient cr'^^ defined by 



l+C„X<Oorl-xsl^<0, 
eif _ (1 +ns)x 



(46) 



(47) 



i^+^nX) 

will become negative and LBS will be unstable. This phenomenon is different from that in Eq. (44). Further, by 
o" = (1/S - l/2)6t and dt = 1, we have 



c? = 



(c?-c?(cr+l/2)x) clU-x/s) 



(48) 

i+inx i-;rs/4 

If S = 2, the correct effective sound speed will be recovered. As indicated in Sec.l, LBS solution in absorbing layers 
is not required to be physical, so it is not necessary for s to be equal to 2. Clearly, because of the collision term in 
LBS, the classical absorbing layer theory can not be used completely and directly. 

As indicated in Sec. 2.2.1, for the type I absorbing term, because of the definition of the effective relaxation s', 
if the condition (23) is satisfied, the instability will grow very fast. In practical applications dealing with acoustic 
problems, because of the nearly vanishing viscosity, s is close to 2 and it is impossible to choose suitable and positive 
X- Therefore, the type I absorbing term is not well suited for absorbing acoustic waves. 

In the same way, we can get the following equation for the type III absorbing term 



d,p' + dxu'j^ 



dtu\ + 



1 + (nX 



dxP 



(1 +ns)x , 

' 1 + ^nX ^ 

(1 +ns)x , 



(49) 



If 1 + („x = 1, in order to keep the sound speed c^ uniform and;^' different from 0, the choice of s must satisfy the 
following equality 

.= ?^. (50) 



However, this requirement is not necessary for the sound speed. The critical value of ;t' for stability is given by 

2 



(1 - 2n + ns)6t 
The positivity constraint on;t' requires the following inequality should be satisfied 

1 



n < 



2-s 



(51) 



(52) 



From the above inequality, it is easy to see that if s is close to 2, the choice of n will become very optional. If we set 
n - 1 12, the constraint on s is natural (s > 0). 

The numerical investigations will be presented in Sec 3.3, which are mainly focused on the type II and type III 
absorbing terms. 

3.2. The linearized LBS with the absorbing terms 

In order to give the linearized LBM with the absorbing terms, we consider the following general form 



fi{x + ViSt, t + 6t)^ fiix, t) + ~s [f^^ip'ix, t), u*(x, t), t) - fi(x, O) + (JfFp' {f/(x, t), u^(x, 0,p*(x, 0, u*(x, t), t)) , (53) 
where the absorbing term is given by 

Ff\pf(x, t), uf (X, t),p% u*(x, t), t)^x [fl"°'\p^i^, t), uf{x, t), t) - /7^V*(x, 0, u*(x, t), t)) . (54) 

Regarding />' (x, f) and u^(x, t) as the uniform reference states, the reference f is defined by 

ff^^f^'^(pf,y^f). (55) 

The fluctuating quantity 6fi is defined by 



Sfi-fi-fl 



(ref) 



(56) 



f(ref) 



f(eq)/ 



f(eq)/ 



Using the definitions of f- and 5fi, and considering /. (p*,u*,f) = / {[fk]Q<k<N)y 'he linearized version of Eq. 
(53) is depicted by 



5fi{x+Vi6tJ+6t) - 6fi(x,t) + s 



r(eq) 



dfr'imo^k^N) 



dfj 



-6: J 

fj=fj I 



5fj+6t5Ff°'^ [pf(x, t), uf(x, t),p*(x, t), u*(x, t), t)) , (57) 



where dPf^ {pf(x, t), u^(x, 0,p*(x, t), u*(x, t), t)) = dPf''^ ({//)o<*<^, {/i)o<i<iv) is calculated by 



SFf''^{{flh<k<NAfk}o<t<N)--X 



dfj 



■ Sfj- 



(58) 



fj=f 



Then, considering a plane wave solution of the linearized equation 



6fj = A jexp[i(k -x-ajt)], (59) 

we get the following eigenvalue problem for the L-MRT-LBM in the frequency-wave number space [4, 5, 8] 



e~"^^f ^M-6i. 



The matrix M in Eq. (60) is defined by 



Mij = dij + s ■ 



dfj 



Jj=fj 



■6,j 



6tx 



^frwfiuj 



dfj 



(60) 



(61) 



fj=fl 



For the sake of comparison with the linearized isentropic Navier-Stokes equations, it is recalled that the analytical 
acoustic modes oj'^ik) and shear modes w*(k) are given by [20] 



Re[w^(k)] = \k\(+cs + |u|cos(k ■ u)). 



lm[w*(k)] 



1 2d-2 



V + rj 



(62) 



2\_t/ 
Re[w"(k)] = |k||u|cos(k-u), 
Im[w^(k)] = -|kpy, 

where v is the shear viscosity and rj is the bulk viscosity. 

3.3. Von Neumann Analysis for D2Q9 model in spectral spaces 

In the above section, the linearized LBS with the absorbing terms was given. In this section, combining the 
analysis presented in Sec. 3.1 and the linearized form (60) in spectral spaces, the numerical dissipative and dispersive 
properties will be analyzed thanks to Von Neumann theory. The parameter 6t for the standard LBS is equal to 1. For 
the sake of convenience, and without loss of gnerality, let us choose s - 1.9 and s - 2m Eqs. (35) and (36). For 
s -2, the molecular viscosity vanishes. One also consider « = m = 1/2. In order to keep the numerical investigation 
identical to the theory given in Sec. 3. 1, let p^ - \ and u^ = 0. The wavenumber in Eq. (59) is defined by 



k ■ cos(0), ;t,, = k ■ sin(0). 



(63) 



Now, let us discuss the stability of the type I absorbing term. In Fig. 2, the dissipation relations for the type I 
absorbing term are shown for two different values of the parameter ;^f. As indicated by the inequality (23), in Fig. 2(a), 
the effective relaxation parameter s' is equal to 2 and the LBS is stable. When the effective relaxation parameter s' is 
larger than 2, it can be observed in Fig. 2(b) that there exist some modes which become unstable. Therefore, type I 
absorbing term is useless for LBS. 

We now investigate the second and third type absorbing terms. First, let us assess the critical value of the absorbing 
strength ;i^ in Eq. (46). If x - 4/?, the effective absorbing strength cr in Eq. (47) will become singular and the 
corresponding dissipation relations are illustrated in Fig. 3. It is easy to observe that under this condition, for two 
values of s {s - 1.99 or ^ = 1.99999), the corresponding LBS is stable. Further, if^ is larger than Ajs, we can observe 
from Fig. 4 that most modes are unstable. These observations validate the theoretical analysis presented in Sec. 3.1. 
The dispersion relations for smaller ;i^ are displayed in Fig. 5. From these results, it is seen that all modes are stable. 
From Fig. 6, the dispersion errors disappear. Certainly, if ;if becomes small, the dispersion error will appears again. 
The suggestion for the choice of ;t' is that for nearly vanishing viscosity acoustic problems, the magnitude of ;t' is in 
the range [0,2]. Theoretically, there exist similar results for the type II and type II absorbing terms. Further, let us 
investigate the type 11 absorbing term by the von Neumann method for several cases . In Fig. 7 , dissipation is shown 
for two choices of;i'. It is easy to conclude that when;^' is larger than the critical ;t' - 4/.?, LBS is unstable; otherwise, 
LBS is stable. 

We consider m^ ^ and {u{, u{) - (0.1, 0). In Fig. 9, u is parallel to k and the dissipation is given for the type II 
absorbing term. Looking at these results, it is clear that the critical j^^ of m-'^ = is suitable for that of u^ + 0. The same 
phenomenon can be observed. However, this kind of phenomenon can not be observed for the type III absorbing term. 
Corresponding results are displayed in Fig. 10. x - ^/■s is not a critical value for the type III absorbing term. Numer- 
ical validation shows that the critical value cifx is smaller than 3/.s. So, it can be concluded that the type II absorbing 
term is the best suited absorbing term for LBS. Meanwhile, we can obtain that for LBS, when the far field u^ is not 
equal to 0, the absorbing terms should damp all possible physical quantities (d,p,da{pu'^) and dp{pu''^u''„ + 6apP*)j in 
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Figure 2: Stability properties of the type I absorbing term: 



:0,n = m= 1/2. 



the macroscopic systems. This is different from the classical absorbing layer for compressible Navier-Stokes equa- 
tions. 
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Figure 3: Stability properties of the type II absorbing term: d = ,n = m = 1/2. 



4. Numerical assessment 

In this section, the type II absorbing term is assessed by solving numerically some classical acoustic problems. 
All numerical investigations were implemented in PalaBos [23]. 

4.1. 2D acoustic pulse 

We first consider a 2D acoustic pulse source. Assuming that the viscosity eff'ect is negligible on acoustic waves, 
the acoustic pulse problem possesses an analytical solution [1, 22]. The initial profile is given as follows 
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Figure 4: Stability properties of tlie type II absorbing term: 9 = ,n = m = 1/2. 
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Figure 5: Stability properties of the type II absorbing term: 6 = ,n = m = 1/2. 



PO = 1 + p/ 
Mo = f/o 
VO =0 



where p/, e, a, r and f/o are defined by 



p/ = eexp(-Q' ■ r^), e = 10 ■', a = ln(2)//7^, r = J(x - xq)^ + (y - yo)^ 



(64) 



(65) 



The computational domain is a [0, 1] x [0, 1] square. And the characteristic length is equal to 1. The parameter b in 
Eq. (65) is equal to 1/20. The exact solution of p/ (if {xo,yQ) - (0,0)) is given by [22] 



pK 






exp I - — I cos(cjrf )Jo(^77)^d^ 
12 



(66) 
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Figure 6: Stability properties of tlie type II absorbing term: 6 = , n = m = \ /2. 0: the exact dispersion relations. 
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Figure 7: Stability properties of the type III absorbing term: 8 = ,n = m = 1 12. 
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Figure 8: Dispersive properties of tlie type II absorbing term: 6 = ,n = m = 1/2. 0: the exact dispersion relations. 
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Figure 10: Stability properties of the type III absorbing term: 6 = ,n = m = 1/2. {Uj^,Uy) = (0.1,0) 



where 77 = -y/(x - Uofp- +y^ + z} and Jo(-) is the Bessel zeroth-order function of the first kind. The adopted lattice 
resolution is 200^. For the sake of convenience, we define the following rescaled time unit 



l/(2f-c,). 



(67) 



Theoretically, when a sponge layer is enforced around the computational domain, after several time units, the 
wave propagate outside the domain and the fluid flow field should converge to the steady reference flow. The time 
for this phenomena is equal to fout = l/( V2f ■ Cj). The ideal state is that the density fluctuation in the computational 
domain vanishes. In order to address this phenomena, one defines the L^ relative error as 



EL<t,x) 



\ 



i:f--(p(x,,o-Prf(x,-,o)2 



2:!lr(pref(x,,o)2 

For the current investigations, pref = 1 in the initial condition (64). 



(68) 
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Figure 1 1 : The schematic figure of the computational domain with the absorbing layer. 
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The computational domain with the absorbing layer is shown schematically in Fig.ll. The normalized (t(x) 
denoted by the formula (34) will be used. The absorbing strength ;t' is chosen equal to Ajs - s, where e = 0.001 denote 
a small positive constant. The new cr for this problem is defined by 

cr(x)=;rc^(x). (69) 

In order to validate influence of the absorbing layer thickness, the thickness of the absorbing layer is chosen equal 
to m ■ b (where m is a positive integer). The current investigations will be compared with the viscosity absorbing 
strategy, which consists in linearly decreasing the relaxation frequency from the physical value to 1 in the sponge 
region. The Reynolds number of this problem is defined hy Re - l/v - 10^ and the viscosity nearly vanishes. For all 
simulations, the Dirichlet boundary condition is used. 

The initial density profile is displayed in Fig. 12. In Fig. 13, the density profiles obtained using the three different 
strategies are shown at different times for a zero mean flow. One can see that without any absorbing strategies a clear 
wave reflection from the boundaries is present. With the linearly decreasing relaxation parameter, there still exits a 
significant wave reflection. With type II absorbing term, only a very weak wave reflection is observed. Especially, 
when the wave propagates across the absorbing layer, it is hardly possible to observe the wave trace. In Fig. 14, a 
quantitative comparison based on EiriUx) is shown. Some clear differences can be observed and type II term is the 
best one. Meanwhile, it can be concluded that using the viscosity damping strategy, the absorbing performance is 
weak. From the given reference line, it can be concluded that the wave in the proposed absorbing layer exhibits a fast 
decay with a decay exponent about -3.5. 

In order to validate the type II absorbing term in the presence of a uniform base flow, the computational results are 
given in Figs. 15 and 16. It is observed in the absorbing layer, the wave has a fast decay exponent about -3. In Fig. 17, 
the error E[2{t,x) is given. From this figure, it is easy to see that when the thickness of the absorbing layers increases, 
the error Eii(t,x) decreases. Meanwhile, a fast decay exponent close to -3.5 is measured. 

These investigations demonstrate that the optimal absorbing term (24) is very eflicient in suppressing the wave 
reflection from the computational domain boundaries. 

Density 
1,00300 

E 1,00200 
1.00100 

1 ,00000 

Figure 12: The initial profile of the 2D acoustic pulse. 



4.2. Validations for time dependent acoustic line source 

In order to validate the absorbing term II for a time-dependent acoustic source, the case of an acoustic line source 
with and without background flow is chosen. The profile of the acoustic line source is described by 

{p(t) - pa+ p'sinilncjt), 
(70) 
U(X) = (M^, Uy), 

where the frequency uj -2. The corresponding wave length is /I = cjl - 0.288675. The schematic domain with 
the absorbing layer is shown in Fig. 18. For present computations, the domain width is, W - AH and the domain 
height is H - 1, respectively, and the width of the absorbing layer is 0.8. The acoustic line source is fixed at the 
left boundary x - (0, y). Periodic boundary conditions are enforced at the top and bottom boundaries. The outflow 
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(a-1) At the time unit T -t 



(c-1) At the time unit T -t 
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Figure 13: The wave propagation at the different times: (a) LBS without absorbing strategy; (b) LBS coupled with the viscosity damping strategy; 
(c) LBS coupled with the optimal absorbing term (24). At the time, the wave crests propagate toward the boundary. At the second time, the waves 
are reflected from the boundaries. The thickness of the absorbing layer is equal to Ab. u' = (f/o, 0) = (0, 0). 
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Figure 14: The L relative error E^2(t,x) '■ -»-: LBS without any absorbing strategy; - < -: LBS coupled with viscosity damping strategy; -o- : 
LBS coupled with the optimal absorbing term (24). The first veilical line indicates waves propagated to the absorbing layers. The second vertical 
line indicates waves propagated to the boundaries, u^ = ([/o, 0) = (0, 0). 



boundary condition is adopted for the right boundary x 
characteristic time scales for zero mean flows: 



(W,y). For the sake of convenience, we introduce two 



T,^W/c„ Tp^llioj-St), 



(71) 



where 6t, To and Tp are the time step, the time of the wave propagation to the right boundary and the time of the 
wave propagation for one wave length, respectively. With the mean flow along x-axis, T,, is defined by 

r,„ = W/(c, + O, (72) 

where m™ is the mean flow velocity along jc-axis. The number of grid points used for this configuration is 200^. 
The absorbing strength ;t' is chosen equal to Aj s - s, where e - 0.001 . Re - 1/v = 10^. 

In Figs. 19 and 20, the computational results for both zero and non-zero background flow are shown. It is observed 
that using the optimal absorbing strategy (24), the wave in the absorbing layer are almost completely damped after 
a short travel distance. Although the given thickness of the absorbing layer is equal to 0.8, the required minimal 
thickness is smaller than 0.3, i.e. about one wavelength A - cjl - 0.288675. Therefore, for this kind of problems, 
the thickness of the absorbing layer can be set to equal to the characteristic wave length or a little larger than the 
characteristic wave length of the acoustic wave. 

4.3. Validations of a dipole advected by the mean flow 

The validation of the optimal absorbing strategy (24) for the vortical flows is very important for the turbulent 
aeroacoustic problems [16], because of the significance of the practical applications. Here, we consider the case of a 
dipole vortex advected by a uniform flow. The dipole is given by 
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Figure 15: The wave propagation at the different times: (a) LBS without absorbing strategy; (b) LBS coupled with the viscosity damping strategy; 
(c) LBS coupled with the optimal absorbing term (24). At the first time, the wave crests propagate toward the boundary. At the second time, the 
waves are reflected from the boundaries. The thickness of the absorbing layer is equal to 4b. u^ = (Uo,0) = (0.1,0) 
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Figure 16: The L relative error E^2{t,x)- -o-: LBS without any absorbing strategy; - < -: LBS coupled with viscosity damping strategy; -o- 
LBS coupled with the optimal absorbing term (24). u^ = ( f/o ■ 0) = (0. 1 , 0) 
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Figure 17: The L relative error E^2(t,x) with respect to the thickness m ■ bof the absorbing layers: -o-: The thiclcness 2b; - < -: The thickness 
2b; ->- : The thickness 3fo; - A -: The thickness 4b. u->' = (Uo, 0) = (0, 0). 
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Figure 18: The schematic computational domain with the absorbing layer 



X 10" 




Figure 19: Comparisons among three different strategies for zero mean flow: (Red solid line) Reference p'; (Dashed-dot line) p' by LBS with the 
optimal absorbing strategy (24); (- < -) LBS coupled with viscosity damping strategy; (- > -) LBS without any absorbing strategies. The sample 
time t = To + IQ ■ Tp. The vertical line indicates the position of the absorbing layer. 
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Figure 20: Comparisons among tliree different strategies for mean flow (u = (u™, »'") = (0.1, 0)): (Red solid line) Reference p'; (Dashed-dot line) 
p' by LBS with the optimal absorbing strategy (24); (- <1 -) LBS coupled with viscosity damping strategy; (- > -) LBS without any absorbing 
strategies. The sample time / = T„, + 10 ■ Tp. The vertical line indicates the position of the absorbing layer. 
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(73) 



The vorticity distribution of the isolated monopoles has the following Gaussian form 



Wo = We 1 - 
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(74) 



If we set We = 299.5285375226, the initial total kinetic energy of the dipolar flow field 
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(x, 0)rdx = 2. 



Considering the mean flow, for the LBS, the initial field is set as 
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(75) 



where Mmean = 0.1. The computational domain is defined as Q = [-1, 1] x [-1, 1]. The parameters in Eq.(75) 
:are taken equal to rg = 0.1, xi - 0, X2 = 0, yi =0.1 and y2 - -0.1. The Reynolds number is set equal to 
Re - 10^ and the computational grid is equal to 400^. The computational domain with the absorbing layer is similar 
with the computational domain indicated by Fig.l 1. The thickness of the absorbing layer is taken equal to 4ro. The 
initial pressure is obtained by solving the pressure Poisson equation and the initial distribution is initialized by the 
method given in [9]. The far field boundary condition is used. The following normalized averaged enstrophy will be 
investigated 
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|w(x, 01 dx. 



(76) 



By monitoring &{t), the absorbing behavior of the proposed optimal absorbing strategy will be validated for vortex 
flows. In all figures, if considering the physical reference velocity is equal to 1, the following time scale is used for 
the LBS computations 

t — t • Mjjiean • v ' ' / 

The initial density field and vorticity field are given in Fig.21. In Fig. 22, the evolution of the density and vorticity 
fields is shown with at three different rescaled times. At the first time shown, the dipole is reaching the absorbing 
layer. At the second time, it is crossing the absorbing layer, while the third snapshot corresponds to a time at which 
the dipole should have left the full computational domain. It is seen that the outlet boundary condition does not radiate 
any spurious wave. The efficiency of type II absorbing layer is measured monitoring time evolution of the enstrophy 
(see Fig. 23). A very fast decay is observed, with a decay exponent close to -12. 

4.4. Flow past two cylinders 

In order to further illustrate the capability of the optimal absorbing layer, the flow past two cylinders with equal 
diameters is considered. The computational domain is drawn schematically in Fig. 24. The centers of the two 
cylinders are located at (lOZ), \QD - \.5D) and (lOZ), \QD + \.5D) where D = 1 is the diameter of the cylinders. The 
thickness of the absorbing layer is taken equal to 2D. The Reynolds number is Re - U^D/v = 5000 (t/^ denotes the 
field field mean velocity magnitude). The Mach number is Ma = 0.0714286. The mesh resolution is equal to D/100. 
For convenience, the following dimensionless time scale is introduced 
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(a)Initial density fields (b) Initial vorticity field 

Figure 21: Initial density field and vorticity field (vorticity is nomialized by the time step St). 



t^t- Uf. (78) 

The reference field velocity u^ = (f/^, 0) and the reference density p^ - \. The three snapshots of the flow fields 
(density, dimensionless velocity magnitude and dimensionless vorticity) are given in Figs.25~27 for t - 49.5, 57.784 
and 66.067, respectively. The very satisfactory quality of the results show that the proposed type II absorbing layer 
can be used to handl complex vortical flows. 

5. Conclusion 

A new absorbing layer technique for the definition of non-reflecting boundary conditions for Lattice-Boltzmann 
methods has been proposed. Thank to an original theoretical analysis of the associated linearized problem, dispersive 
and dissipative features of the method have been analyzed. An interesting result is that the associated sponge layer 
equation for the associated macroscopic problem is more complicated that sponge layer equations directly defined 
on macroscopic quantities. By recovering the corresponding macroscopic equations, the macroscopic effect of the 
absorbing terms is identified and the critical absorbing strength value is formulated. Performing mathematical analysis 
in the spectral space, the influence of critical parameters are investigated. From this analysis, an optimal absorbing 
strategy is obtained. By simulating two basic acoustic problems, the optimal absorbing strategy is validated and the 
corresponding decay exponent is obtained. Numerical results show that the optimal absorbing strategy is also powerful 
for damping vortices. A very high damping exponent for the enstrophy is observed. The results demonstrate that the 
optimal absorbing strategy Type II is very effective for coping with vortical flows. 
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Figure 22: Evolution of density field and vorticity field (vorticity is normalized by the time step St). 
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Figure 23: Time dependent normalized enstropliy G(t). Tlie vertical lines indicate the time at which the dipole moves across into the absorbing 
layer. The slope of the reference oblique line is equal to -12. 




Figure 24: The schematic computational domain with absorbing layer for the flow past two cylinders. 
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(c)F= 66.067 

Figure 25: The snapshots of the density fields. 




(c)f= 66.067 

Figure 26: The snapshots of the dimensionless velocity magnitude fields 
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(c)f= 66.067 

Figure 27: The snapshots of the dimensionless vorticity fields 



Appendix A. Transformation matrices 

Appendix A.l. For two dimensional MRT-LBM with 9 discrete velocities 

The transformation matrice for MRT-LBM with 9 discrete velocities is described by 



M 
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-1 
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(A.l) 
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Appendix A.2. For three dimensional MRT-LBM with 15 discrete velocities 

The transformation matrice for MRT-LBM with 15 discrete velocities is described by 



M 



1 1 
1 -1 





1 



-1 





1 1 1 

1 -1 
1 1 



1 -1 



-2 -1 -1 -1 -1 -1 -1 

16 -4 -4 -4 -4 -4 -4 

1-1 
-44001 1 

0-4 4 1 1 
-1 -1 -1-10 

1 1-1-10 
1-1 
1 1 
1-1 






-4 


4 























2 


2 






















































1 -1 1 


-1 1 


-1 


1 1 -1 


-1 -1 


-1 










1 -1 1 


-1 1 


-1 


1 -1 1 


1 -1 


-1 


1 1—1 


— 1 —1 


— \ 






1 






1 



-1 -1 -1 -1 
1-1-1 1 






-1 

1 

-1 



1-1-1 1-1 



1 



1 -1 



(A.2) 



Appendix B. The macroscopic equations with the absorbing layer for the BGK model 

Appendix B.l. The classical absorbing layer 

The BGK model with the absorbing layer is given by (0 < /, j < N) 

Mx + Vidt, t + 5t)^ fi(x, t) + s {f'p\p*{x, t), u*(x, t\ t) - //x, t)) + 6tx {ff''\p^(x, t), u-^(x, t), t) - //x, t)), (B.l) 

where p*(x, t) and u*(x, t) are the equilibrium density and velocity, respectively. Eq. (B.l) can be rewritten as follows 

fi(x + ViSt, t + 6t)^ fi(x, t) + s' {fi'"^\p*(x, t), u*(x, t), t) - fiix, t)) + 6tFi(pf(x, t), u^(x, t),p*,u*(x, t), t), (B.2) 

where s' = s + 6tx and F,(u/(x, f), u*(x, f), is defined by 

F,(p-'(x, 0, u-^(x, t\p\x, t), u*(x, 0, t)^x [ff''\pH^, t\ uf{x, t), t) - /;.^^^V*(x, 0, u*(x, 0, t)) . (B.3) 

Eq. (B.2) can be regarded as the new LBS with a special external force term. In order to obtain the macroscopic 
equation corresponding to Eq. (B.2), the classical Chapman-Enskog derivation [17, 21] is used. The following 
expansions are introduced [2 1 ] 

/. . ff> + eX<^' + e'f' + . . . , (B.4) 

d, = ed,, + e%, da = edoa, F, = eF°. (B.5) 

By the multi-scale expansion, we have 



e' : f' ^ ft%\u\x,t),t) 



f(0) rO 



e' : Du/r - F'/ = -^' m 



m 



,2 ^(0) 



f(2) 



e' : d,jr + D,Jl" + St/2Difl"> = -s'/6tf^^>. 



(B.6) 
(B.7) 
(B.8) 
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where Do, - dti + ciada- Eq. (B.20) can be expressed as follows 

Here, we define the macroscopic quantities {p* , pu*^) as follows 

p* =2_^fi + nSt 2_^ Fi, p*M* = 2_) Ciafi + mdt 2_j CiaFi. 

(■ / i i 

It is easy to get the following relation 

2 f' = -nSt X Fl X c^f^ = -m6t ^ c,F°. 



Ak) 



So, f^ has the following properties 



and 



At the scale fi, we have 



Zr = o,^>i, 

i 

5„p* + 5o.(p*<) = X ^'^ + "^' Z ^'°' 



5„ (p*M* ) + dop (;r^y) = ^j '^'■»^" + ""' 2j ^'■'^^'°' 



.(1): 



where the second-order tensor n^' is defined by 



M) 



At the scale t2, we have 



111 *** .*c 

%/3 = P «a« -/* +/? Safl 



6t 



d„p* +6til--\\ -nd„ Yj^i- '^^Oa Yj ^-^f + T K> 2 ^^ + ^''" Z ^-^' 



5f,(p*M*)+ 1 



+ 
2 



i ifi ) V 1 ifi 



•(I) 



Now, let us calculate Yjifi ciaCipf^ . From Eq. (B.7), we have 



St 



2_j ^iaf^'Pfl^^ = ^la 2j ^iaCilsff''^ + doy 2_j CiaCipCjyfl^^ - ^ CiaCi/jF^ . 



(B.9) 
(B.IO) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 
(B.15) 

(B.16) 

(B.17) 
(B.18) 

(B.19) 



Then, we have 



where/ =c^,(l/s' - l/2)6t andcr = (l/s' - 1/2) (Jf. 
Eq. (B.20) can be rewritten as follows 



1 - y 2_i ^io^ipfl^^ ^ vp(dQaUp + dapuj - a-diy,(pulupU^) - cr^j CiaCipF^i- (B.20) 



d,^(p*ul) - dopiv'pidoaUp + dopuj) + crdopdoyip* ulupU^) + crdop ^ CiaCipF^- + 



1- 



-m6td,„ Y 



c F' 



^ 6t^ 

^2 



ifi 
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= 0. 



(B.21) 



Combining Eqs. (B.14), (B.15), (B.17) and (B.21), ignoring the term of 0(Mfl^),we obtain 



dtp* + d„(p*ul) = /] ■P'l + «■?' 2_1 ■^' + ^M 1 ~ 7 "'5' Zj ■^' "*" '"^" Zj '^''"■^'' ~ T r ' Zj ■^' "*" '^" Zj '^'"■■^' 



, (B.22) 






6t\l-- 



md, 2_] CiaFi 



^ 6t^ 



(B.23) 



From the defination of F,, we have 

2j ^'' = AfOo-^ - P*). 2j '^''^^' ^ ^ (^"^"^ ~ '°*"«) ' Zj ^i"^>P^> ^ ^ {p^uiu^ + pfdap - p*ulup - p*6ap) . (B.24) 

/ / i 

So, by Eq. (B.23), we have 

d,p*+da(p*ul) = {\+ns')x(pf-p*)+(n5til - '-^-^Axd.ipf -p*)+x(rnSti\ - ^j - |J5„ (p-^«{ -p*<) , (B.25) 

dt(p*ul) - dp(v'p*(dlup + dpuj) + dp [p*ulup + p*5apj = (1 + ms')x [p^ui - p*u^ + 
\m6t y-^]-^ \xdt [p^ui - P*u*a) - (o" + y JAf^/^ (p-^MffW^ + P^^ap ' p'K^p ' P*6ap) ■ (B.26) 



2/ 2 
Let (,„ and ^m defined by 



s'\ 5t s 



6t 



Eqs. (B.25) and (B.26) are rewritten as 



d,p' + daip'ul) = (1 + ns')x(p^ - p") + ^nXdtip^ - P") + (mX^a {p^ ui - p* U^) , 

d,(p*ul) - dp{v'p*{daUp + dpuj) + dp {p*ulu*p + p*6ap) = (1 + ms')x {p^ul - p*u^ + 

(wXdt [P^ui - P*Ma) - (o" + -l^jXd/S [p^uUp + P^Sap - p*ulUp - p*6ap) ■ 



(B.27) 



(B.28) 



(B.29) 



Appendix B.2. The absorbing layer based on the equilibrium distribution functions for the BGK model 
The BGK model with new damping terms is given by (0 < /, j < N) 

fix + Vi6t, t + 6t)^ fix, t) + s [ff^p*, u*(x, 0, - /,(x, t)) + StFi (uy(x, t), u*(x, t), t)) , (B.30) 

where u*(x, t) is the equihbrium velocity. F (uy (x, t), u*(x, t), t) is defined by 

Ff'{pf,M\x,i),p\M\x,i\t) ^x[ff''\p*\^\^^\t)-ff°'\p\^y\^^),t)). (b.31) 

From Eq. (B.30), the obtained macro equations are similary to Eqs. (B.28)-(B.29). The parameters v, („ and ^,„ are 
defined as follows 



S - S, V - c. 



(^-^)^r,^„ = «c5r(l-^)-f,^„ = ,«.r(l-^) 



5\ 6t 

2' 



(B.32) 
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Appendix B.3. The absorbing layer based on linear damping terms 

In Appendix B.l, the absorbing terms are specified based on the equilibrium distribution functions. From Eq. 
(B.29), it is known the obtained damping terms in the recovered Navier-Stokes equations involves a nonlinear damping 
term (the second order terms of the velocity u' and u*) in momentum equations. The following model is propsed for 
elinimating this term 

fiix + vidt, t + St)^ fi(x, t) + s (/f ^V*, u*(x, t), t) - fj(x, O) + 6tFf' ^^> (u^(x, t), u*(x, t), t)) , (B.33) 

where u*(x, t) is the equilibrium velocity. F ' '^'' (u/(x, t), u*(x, t), i) is defined by 

Ff' '''\pf, uf(x, t),p% u*(x, t), t)^x [t' ''V^, uf(x, t), - /^- "''>(p\n*(x, t), t)) . (B.34) 

Then, along the routine in Appendix B.l, the following macroscopic equations are obtained 

d,p* + da(p*u*^) = (1 + ns')x(p^' - p*) + CnXdtip^ - P*) + Uxda {p^ui - p*ul) , 
dtip*K) - dfj(vp*(d„u*^ + dpuD) + dp {p*uy*^ + p*Sap) = (1 + ms')x [p^ui - p*ul) + ^^xd, [p^ui - p*ul) , 

where the parameters v, („ and („, are defined as follows 



(B.35) 



S, V = c. 



(l-i).r,^„.«c5r(l-^)-f,^„ = ,«.r(l-^)-f. (B.36) 



References 



[1] H. Xu, P. Sagaut, Optimal low-dispersion low-dissipation LBM schemes for computational aeroacoustics, J. Comput. Phys 230 (13): 5353- 

5382(2011) 
[2] F. Dubois, P. Lallemand, Towards higher order lattice Boltzmann schemes, J. Stat. Mech. Theory E, (2009) P0600.6 
[3] J. M. Buick, C. A. Created. D. M. Cmpbell, Lattice BGK simulation of sound waves, Eurohys. Lett. 43 (2) (1998) 235-240. 
[4] S. Marie, D. Ricot, P. Sagaut, Comparison between lattice Boltzmann method and Navier-Stokes high order schemes for computational 

aeroacoustics, J. Comput. Phys. 228 (2009) 1056-1070. 
[5] D. Ricot, S. Marie, P. Sagaut, C. Bailly, Lattice Boltzmann method with selective viscosity fiher, J. Comput. Phys. 228 (2009) 4478-4490. 
[6] J. M. Buick, C. Spectra and L. Buckley, C. A. Created, Lattice Boltzmann BGK-simulation of non-linear sound waves: The development of a 

shock front, J. Phys. A: Math. Cen. 33 (2000) 3917-3928. 
[7] S. Chen, C. Doolen, Lattice Boltzmann method for fluid flows, Annu. Rev. Fluid Mech. 161 (1998) 329. 
[8] P. Lallemand, L. S. Luo, Theory of the lattice Boltzmann method: Dispersion, dissipation, isotropy, Galilean invariance, and stability, Phys. 

Rev. E. 61(6) (2000) 6546-6562. 
[9] H. Xu, H. B. Luan, Y. L. He, W. Q. Tao, A lifting relation from macroscopic variables to mesoscopic variables in lattice Boltzmann method: 
Derivation, numerical assessments and coupling computations validation. Computers & Fluids, doi:10.1016/j.compfluid.2011.10.007, 2011. 
[10] H. Xu, W. Q. Tao, Y. Zhang, Lattice Boltzmann model for three-dimensional decaying homogeneous isotropic turbulence. Physics Letters A, 

373 (15) (2009) 1368-1373. 
[11] D. d'Humieres, I. Ginzburg, M. Krafczyexact solution of k, P. Lallemand, L. S. Luo, Multiple-relaxation-time lattice Boltzmann models in 

three dimensions, Phil. Trans. R. Soc. Lond. A, 360 (2002) 437-451 
[12] T Colonius, Modeling artificial boundary conditions for compressible flow, Annu. Rev. Fluid Mech. 36 (2004): 315-345. 
[13] D. J. Bodony, Analysis of sponge zones for computational fluid mechanics, J. Comput. Phys. 212 (2006) 681-702. 
[14] M. Israeli, S. A. Orszag, Approximation of radiation boundary conditions, J. Comput. Phys. 41 (1981) 115-135. 
[15] J. B. Freund, A simple method for computing far-field sound in aeroacoustic computations, J. Comput. Phys. 157 (2000) 796-800. 
[16] C. A. Wagner, T. Hiittl, P. Sagaut, Large eddy simulation for acoustics. New York: Cambridge University ftess, 2007. 

[17] Z.L. Guo, C.G. Zheng, B.C. Shi, Discrete lattice effects on the forcing term in the lattice Boltzmann method, Phys. Rev. E. 65 (2002): 046308. 
[18] D.R. Durran, Numerical methods for fluids dynamics with applications to Geophysics, Springer New York Dordrecht Heidelberg London, 

2010. 
[19] R.R. Renaut, J. Frohlich, A pseudospectral Chebychev method for the 2D wave equation with domain stretching and absorbing boundary 

conditions, J. Comput. Phys. 124 (1996) 324-336. 
[20] L.D Landau, E.M. Lifshitz, Fluid Mechanics, second ed., Oxford: Pergamon, 1987. 
[21] J.M. Buick, C.A. Created, Gravity in a lattice Boltzmann model, Phys. Rev. E, 61 (5) (2000): 5307-5320 
[22] C.K.W. Tam , J.C. Webb, Dispersion-relation-preserving finite difl'erence schemes for computational acoustics, J. Comput. Phys. 107 (1993) 

262-281 
[23] User's Guide, http://www.lbmethod.org/palabos/. 

32 



